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We consider type I+II seesaw mechanism, where the exchanges of both right-handed neutrinos 
and isotriplet Higgs bosons contribute to the neutrino mass. Working in the left-right symmetric 
framework and assuming the mass matrix of light neutrinos and the Dirac-type Yukawa couplings 
to be known, we find the triplet Yukawa coupling matrix /, which carries the information about the 
masses and mixing of the right-handed neutrinos. We show that in this case there exists a duality: 
for any solution /, there is a dual solution / = m^/vL — f, where vl is the VEV of the triplet Higgs. 
Thus, unlike in pure type I (II) seesaw, there is no unique allowed structure for the matrix /. For n 
lepton generations the number of solutions is 2". We develop an exact analytic method of solving 
the seesaw non-linear matrix equation for /. 
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PACS numbers: 14.60.Pq, 14.60.St, ll.30.Er 

1. Introduction. The discovery of neutrino oscil- 
lations in atmospheric, solar, reactor and long-baseline 
neutrino experiments [l| gave an unambiguous evidence 
for neutrino mass and mixing and spurred the activity 
in neutrino mass model building. Among the possible 
mechanisms of neutrino mass generation, arguably the 
simplest and most attractive one is the seesaw mecha- 
nism 0, Q , which explains the smallness of the neutrino 
masses through the exchange of superheavy particles. In 
addition, it has an elegant built-in mechanism of gener- 
ation of the baryon asymmetry of the universe ~ baryo- 
genesis via leptogenesis 

In its simplest version, the seesaw mechanism extends 
the standard electroweak model by assuming the exis- 
tence of S'C/(2)i-singlet right-handed (RH) neutrinos Nji 
with a bare Major ana mass matrix Mji. The RH neu- 
trino masses, not being protected by the electroweak 
symmetry, are naturally large ~ at Planck or GUT scale. 
The masses of light neutrinos, produced in this (so-called 
type I) seesaw mechanism, are 



(1) 



where mo = vy, v ~ 174 GeV being the electroweak 
VEV and y the Dirac-type Yukawa matrix describing 
the coupling of lepton doublets to Nj^ and the standard 
Higgs doublet. Alternatively, the masses of light neutrino 
can be generated through their interaction with a heavy 
S'?7(2)i-triplet Higgs scalar A^, which can acquire an in- 
duced VEV Vl through its coupHng to the Higgs doublet. 
This (so-called type II) seesaw mechanism Q yields 
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where is the Yukawa coupling between the lepton dou- 
blets and Al, and Ma is the mass of the isotriplet Higgs. 

The RH neutrinos Nr, being electroweak singlets, 
are actually aliens to the standard electroweak model. 



though can be freely added to it. They are much more 
natural in models with left-right (LR) symmetry or 
in SO{10) GUT 0], which can be broken down to the 
standard model through the LR-symmetric stage. The 
LR symmetry provides a very natural explanation of the 
observed maximal parity and C violation in low-energy 
weak interactions and is therefore likely to be present at 
some level in the final theory. 

In this class of models, the RH neutrinos are in dou- 
blets In of SU{2)ji, and their mass is generated by the 
VEV VR of the S't/(2)fl-triplet Higgs field A^, through 
the Majorana-like Yukawa coupling. LR symmetry then 
implies that there should also exist an S'C/(2)L-triplet 
Higgs field A^, and the Yukawa couplings of the triplet 
Higgs scalars to leptons have the form 

h (lie IT2 Al II) + /fl (^f C iT2 Ar II) + h.c. , (3) 
T 

where 11 = CIr, , Al^r = tA.l^r and C is the charge 
conjugation matrix. The Dirac-type Yukawa matrix y 
comes from the coupling of the lepton doublets / l and Ir 
with the bi-doublet Higgs field In this class of models, 
light neutrino masses naturally get both type I and type 
II contributions, with Mr = vr/r. 

In addition to the gauge LR symmetry, the models of 
this kind are usually assumed to possess a discrete LR 
symmetry, which is broken at a scale vlr that may or 
may not coincide with the SU{2)r breaking scale vr. 
There are essentially two ways in which the discrete LR 
symmetry can be introduced. In the first one, the left- 
handed and RH fermions are interchanged, while the 
Higgs fields undergo A^ ^ Ap, $ ^ $^ This yields 



Jl - f*R: 



(4) 



where the last equality holds when the VEVs of the 
bi-doublet $ are both real. The second realization re- 
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quires the invariance w.r.t. Il ^ 11, A/, <-^- = A^^, 
$ ^ = (j)^^ i.e. is essentially the charge conjugation 
symmetry. It leads to 



yields 



y^y 



(5) 



While both implementations of the discrete LR symme- 
try are possible, the second one is more natural in the 
context of 50(10) unified theories, where it is an auto- 
matic gauge symmetry 7] . In what follows we will adhere 
to this latter realization, i.e. will be assuming the equal- 
ities in eq. © to hold. The light neutrino mass matrix 
in this LR-symmetric seesaw mechanism is then given by 



nil. 



(6) 



For the realization Q of discrete LR symmetry, type I 
term in eq. © would contain (/*)~^ rather than /^^. 

2. Seesaw and duality. We will pursue a bottom- 
up approach: assuming rrii, and y to be known, we solve 
eq. lO for the n x n matrix /, where n is the number of 
lepton generations. Note that the matrices and / are 
in general complex symmetric, whereas y is symmetric 
due to the assumed discrete LR symmetry 

An examination of eq. lO reveals the existence of its 
following duality property: if / is a solution, so is 



/ 



nil, 

VL 



(7) 



provided that the matrix y is invertible. This can be read- 
ily verified by the direct substitution of |7J| into eq. ©. 
Notice that exactly the same duality holds in the case of 
the realization Q of the discrete left-right symmetry. 

The duality implies that, given and y, there is no 
unique solution for the triplet Yukawa coupling /. As we 
shall see, for n fcrmion generations there are 2" matrices 
/ that, for a given y, result in exactly the same mass 
matrix of light neutrinos ni^. It should be stressed that 
this duality is a unique property of the LR-symmetric 
type I-I-II seesaw: no such duality exists in the cases of 
pure type I or II seesaw, nor in the case when the two 
seesaw contributions are unrelated. 

In LR symmetric models with non-minimal particle 
content, the presence of extra particles may complicate 
the structure of the seesaw. In that case, in general there 
is no duality (at least in the form discussed in this paper). 

We shall now consider the illustrative examples of one 
and two lepton generations as well as the realistic three 
generations case. 

3. One lepton generation. In this case m^, y and / are 
merely complex numbers. Eq. JBJ is a quadratic equation 
for / with solutions 



f± 
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The quantities f± satisfy /- +/+ = m^/vL, i.e. are duals 
of each other. In the limit jm^p ^ 4|y^w^Wi/uii;|, eq. © 



y^v 
m^VR 



U^^+y^. (9) 
VL m„VR 



This limit leads to the dominance of type I or II see- 
saw in the light neutrino mass, depending on which so- 
lution, /_ or f^, is realized in nature. For ^ 1"^!^!, 
the solution violates the perturbative unitarity and 
must be discarded, so that only /_ can be acceptable; a 
necessary condition for either solution to be physical is 
\vr\ > \v'^y'^/nii,\. 

When the two terms under the square root in eq. © 
are of the same order, the contributions of type I and 
type II terms to are also of the same order (mixed 
seesaw). Finally, in the limit \ni^\'^ ^ Mv'^^'^vl/vrI, one 
has ml ~ —nil^ and |to^'^^| 3> \ni^\, i.e. in this case 
there is a strong cancellation of the two terms in eq. 

Consider now n fermion generations with no mixing 
(which means that / and y are diagonal in the same 
basis). This case can be described as n replicas of the 
one-generation case: there exist n pairs of solutions f±^^ 
(i — 1, ...,n) to n quadratic equations, and each solution 
for the matrix / corresponds to picking one value f^^ 

or from each of n pairs. The number of solutions is 
thus 2" . Switching on the mixing will change the nature 
of the solutions (which will no longer be diagonal), but 
not their number 

4- Two lepton generations. In this case m^, y and / 
are symmetric 2x2 complex matrices. For definiteness, 
we shall consider the (i'^, Vr) sector of m^,. The matrices 
and y are given by 



1 

F 



fss —.f23 
— f23 /22 



y 



2/r2 2/r3 



(10) 



where F = detf and j/t-2 = y^is- Defining x = vlVr/v'^, 
one can rewrite eq. in components as 

/33y^2 - 2/23y/x22/M3 + /222/^3 > 



/33y^3 ~ 2/23j/^3yr3 + /22y?3 ' 



xF{f22 - TO^^) 

xF{f33 ~ mrr) - JSiy^3 

xF{f23-mf,r) = /33yM2yM3-/23(2;M2yr3 + y^3) 

+ /22yM32/r3 , (H) 



where m = m^jvL. This is a system of coupled cubic 
equations for the matrix elements of /. Usually, solv- 
ing systems of this kind presents serious difficulties; one 
has to resort to numerical methods, whose application is 
complicated by the fact that non-linear equations have 
multiple solutions, and their number is often not known 
in advance. We shall, however, show now that the system 
of equations Hll|) admits a simple exact analytic solution. 
To find this solution, we develop here a method that may 
also prove to be useful for solving similar systems of cou- 
pled non-linear equations appearing in different contexts. 
Let us rescale all the matrices of interest according to 



^J - , y^J 



^yi 



(12) 
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where A is an as yet arbitrary complex number. The 
scahng law was chosen in such a way that in terms of the 
primed variables the system of equations for //^ has the 
same form as eq. (|llf) . Next, we fix the value of A by re- 
quiring F' = detf — 1. The system of equations for f^j 
then becomes linear and can be readily solved. Express- 
ing the primed variables back through the unprimed ones 
and substituting them into the condition F' = 1, one ob- 
tains a 4th order characteristic equation for A. Solving it 
completes the solution of the problem. 

Since the matrix y is symmetric, one can go, without 
a loss of generality, to the basis where y is diagonal: y — 
diag{y2, ys)- The solution for the matrix / then takes 
the form 

J f xXnifj^fj, + ylrrirT m^.rixX - y2y:i)\ 

{xXy - ylyl \ ■■■ xXmrr + vlm^.^,) ' 

(13) 

where A has to be found from the characteristic equation 

[{xXf - ylyl] -X [det m{xX - y2yifxX 

+ {mt,t,yi + mrrV2f{xXf]={). (14) 

Eq. (|14|l is quartic in A and so has in general four complex 
solutions, leading to four solutions for the matrix /, as 
expected. 

Rewriting eq. © as x/ = —yf ^^y and taking the de- 
terminants of both sides, one finds x^FF — j/|y|, where 
F = det /. Since the condition F' = 1 corresponds to 
F = A, this equation leads to x'^XX = With the 

help of this relation, it is straightforward to check that 
the four structures of / defined by eqs. H13|) and H14|) 
form two dual pairs, and it is also easy to express the 
four solutions of eq. (|14|l in a closed form: 



xXi ~ — [x det m-f'r±±-\/2(dct m)'^x'^ + Akx + 2xr± det m 

(15) 

where 

r± = ±\J (det m)2x^ -I- Akx + IGy^y^ , 

k = ml^yl + 2ml^y2y3 + ml^yl ■ (16) 

Eq. (|15|l gives two pairs of dual solutions, one of them 
corresponding to and another to r_ . In each dual pair 
the solutions differ by the sign in front of the radical. 

From eq. H13|) one can see that when a solution Ai 
satisfies |a;Ai| ^ \yiyj\ (hj — 2,3), one obtains /i ~ m, 
which corresponds to type II seesaw case. Then the dual 
solution xX2 = xXi = y^yi/ (xXi) has modulus <C lyiyjl 
and the corresponding matrix /2 = /i takes the form 
obtained in type I seesaw case. Notice that eqs. (I15II and 
(|16|1 imply that the condition |a;Ai| ^ \yiyj\ can only be 
satisfied in the limit Imafsm^sl ^ Mytyj/^l {^^tPtI^^ = 
/J,, r). This condition ensures the existence of solutions 
with the dominance of one seesaw type, in analogy with 
the one generation case. However, in general not all four 
solutions correspond to the one seesaw type dominance in 



this limit: if | detm| ^ A\yiyj/x\, the remaining two dual 
solutions /a and are of mixed type. When \mai3'm-ys\ ^ 
A\yiyj/x\, in general all four solutions are of mixed type. 
A similar classification applies to the 3-generation case. 

5. Three lepton generations. In this case /, y and 
m — m^/vL are complex symmetric 3x3 matrices. As 
in the 2-generation case, we go to the basis where y is 
diagonal: y = diag{yi, j/2, ya)- Quite analogously to the 
derivation of eq. Hll|l , one can obtain from equations 
for the matrix elements of / and its dual /: 

xPifij - iTT'tj) = ytyjFtj , (17) 
xFifij ~ n^tj) = -xFfij = yiyjFi-i , (18) 

where F = det f,F = det / as before, and 

Taking the determinants of both sides of the equation 
xf ~ ^yf~^y, one obtains 

x^FF = -ylylyl . (19) 

Eq. H17|l is a system of six coupled quartic equations 
for the elements fij (note that F is cubic in /y in 
the 3-generation case). Since the right-hand sides of 
these equations are quadratic rather than linear in /y , 
a simple rescaling would not immediately linearize the 
system. However, using the dual system of equations 
(|18|l one can write for these right-hand sides yiyjFij — 
—xFfij + yiyj(Tij — Mij). The resulting system of equa- 
tions can now be linearized by a rescaling, similar to 
that in eq. H12|) . except that the scaling factor is now 
A^/"^ rather than \/A. Once again, we fix A by requiring 
F'(A) = 1; eq. ^ then yields F' = - {y'iy2y3f / x^ . The 
linearized system is 

- {y'MfUl, ~ = ^yly-m, - ^4) • (20) 

The characteristic equation for A is now of 8th order and 
its solutions lead to four pairs of dual solutions / and / 
for the system (|17|l . This 8th order equation and the re- 
sulting general analytic solution for / are rather lengthy 
and will be given elsewhere; here we present the simpli- 
fied case yi 0. 

The case \yi\ <C |?;2,3| is actually physically well mo- 
tivated, in view of the smallness of the masses of first- 
generation charged fermions. Notice that for yi — the 
matrix y is not invertible, so that the proof of the duality 
of solutions for / breaks down. However, one can check 
that the solutions of eq. H17|l for yi = are the same as 
the ones obtained by setting ?/i = in eq. (|2()|l . In other 
words, one can safely use the duality and take the limit 
j/i ^ in the final result. The solution for the symmetric 
matrix / is then given in components as 

/ll = , /l2 = mefj, , /13 = Trier , 

/23 = [rn,r + ) /d2 , 
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/22 — 

/33 = 




where 



di = 1 - 



d2 = 1 



X 



X 



(22) 



In the considered hmit j/i ^ the characteristic equa- 
tion for A reduces to 

A'* {[(sA)^ - ml^ylulf - x [detm(.xA - meey2V3)^xX 

+ (Af222/2 + A/332/3)' (xA)2] } = . (23) 

A comment on the duahty and multipUcity of solutions 
in the 3-generation case is in order. If a solution A of the 
general characteristic equation is non-zero in the limit 
yi — > 0, then A = — yiylj/i/l^'^'^) ^ Oj which means that 
the determinant of / vanishes. Thus, the dual of any so- 
lution that is finite for yi — *■ becomes singular and must 
be discarded. Therefore, in this limit there are only four 
(rather than eight) solutions with no duals. The corre- 
sponding values of A are the zeros of the factor in curly 
brackets in eq. H23|l . which is quartic in A. From the com- 
parison with eq. I|14|l a strong connection with the pure 
2-generation case is evident, so that a different duality 
among the four remaining solutions exists: if A satisfies 
eq. (f^ . so does A = 2/|2/|mgg/(a;^A), and it corresponds 
to / = TO — /, where rhafs = map -\- meamc-p/mee- There 
are two pairs of such solutions. 

6. Concluding remarks. We have revealed and an- 
alyzed an interesting duality property of the seesaw 
mechanism with contributions of both RH neutrinos and 
isotriplet Higgs bosons in the presence of discrete left- 
right symmetry. In particular, it has been shown that 
for a given mass matrix of light neutrinos TOjy and Dirac- 
type Yukawa coupling matrix y, there are multiple solu- 
tions for the Majorana-type Yukawa coupling matrix / 
and thus for the matrix of RH neutrinos Mr. The num- 
ber of solutions, however, does not exceed eight (for three 



lepton generations), and with the help of the developed 
here formalism they can be readily analyzed one by one. 
In contrast to this, in models where type I and type II 
contributions to neutrino mass are unrelated, there are 
in general infinitely many possible decompositions of to^ 
into these two parts. 

The discrete LR symmetry of the underlying theory 
must be broken at some scale vlr, and the renormaliza- 
tion group evolution below this scale may result in a vio- 
lation of conditions lO at lower energies. The corrections 
to the matrix elements of y and / depend logarithmically 
on the ratios of the masses of RH neutrinos and Higgs 
triplets and are suppressed by loop factors and possibly 
by small couplings. We checked numerically that the re- 
construction of the matrix / as performed in this paper 
remains accurate at a percent level if the LR-violating 
corrections to the matrix elements are of the order of 
percent. This stability, however, may be lost if small 
matrix elements receive corrections proportional to the 
large ones. This issue requires a dedicated study, which 
is beyond the scope of the present paper. 

Possible applications of our results include the bottom- 
up reconstruction of the mass matrix of heavy RH neu- 
trinos and model building. The bottom-up approach re- 
quires the knowledge of the light neutrino mass matrix 
and of the Dirac-type Yukawa couplings y. The former 
is partially known from the experiment (absolute mass 
scale, 1-3 mixing angle and CP violating phases have not 
been measured yet), whereas for y one needs both data 
and additional theoretical assumptions, such as quark- 
lepton symmetry or Grand Unification. Model builders 
can use our results to immediately find out the allowed 
structures of the triplet Yukawa couplings / that can suc- 
cessfully reproduce the low energy phenomenology, which 
may help them to systematically look for the underlying 
symmetries of the mechanism of lepton mass generation. 
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